THE DEFECT OF GENERALIZED FOURIER MATRICES 



TEODOR BANICA 



Abstract. The N x N complex Hadamard matrices form a real algebraic manifold Cn- 
We have Cn = Mn{T) \/NUn, and following Tadej and Zyczkowski we investigate 
here the computation of the enveloping tangent space ThCn = ThMnCT) n ThVNUn, 
and notably of its dimension d{H) = dim(Tf/CAr), called undephased defect of H. Our 
main result is an explicit formula for the defect of the Fourier matrix Fq associated 
to an arbitrary finite abelian group G = Z^i x . . . x Z^r^ . We also comment on the 
general question "does the associated quantum permutation group see the defect" , with 
a probabilistic speculation involving Diaconis-Shashahani type variables. 
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Introduction 

A complex Hadamard matrix is a square matrix M G Miv(C), all whose entries are on 
the unit circle, \Hij\ = 1, and whose rows are pairwise orthogonal. If we denote by T the 
unit circle, and by Un the unitary group, the set formed by such matrices is: 

Cn = M^(T) n VnUn 

Observe that Cat is a real algebraic manifold. In general Cat is not smooth, nor it is a 
complex algebraic manifold. Its structure is in fact very complicated. In order to describe 
it, best is to introduce the quotient set Cn — ^ consisting of complex Hadamard 
matrices under the standard equivalence relation between such matrices, obtained by 
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permuting rows and columns, or multiplying them by complex numbers of modulus 1. At 
N = 2,3,4,5 this set is well understood, thanks to the results of Haagerup in [16]. 
At = 6, however, this set En is only known to contain: 

(1) Two tori T^, intersecting at the Fourier matrix Fq. 

(2) A circle T, coming from the Haagerup matrix Hq. 

(3) An isolated matrix Tg, discovered by Tao in [29] . 

(4) A circulant matrix Cg, found by Bjorck and Froberg in [12j . 

(5) An algebraic curve Sq, found by Beauchamp and Nicoara in [^j. 

(6) Some other components, all quite poorly understood, see [25]) ^li- 
lt is quite unclear, starting from this data, on how to develop a systematic study to 

Cat. An interesting idea, however, comes from the paper of Tadej and Zyczkowski [28j . 
Motivated by some early work of Nicoara p3] and by questions regarding unistochastic 
matrices |10], they computed the "enveloping tangent space" at H E Cn, given by: 

ThCn = ThMn{T) n ThVnUn 

The description of ThCn found in [28] is of course quite a theoretical one, in terms of 
solutions of a certain system of linear equations associated to H. Of particular interest is 
the dimension of the space of solutions, called "undephased defect" of H: 



d{H) = dimiTnC^ 



N) 



Tadej and Zyczkowski did as well in [28] a number of defect computations. Their main 
result there is that, for the Fourier matrix of order A^ = . . .pg% we have: 



d{FN)=Nl[(l 

1=1 ^ 



tti 

Pi 



Back to the general case now, knowing ThCn is of course quite far from understanding 
the local structure of Cat around a given point H E Cn- As pointed out by Barros e Sa 
and Bengtsson in [8], even in the case of the Fourier matrix F12 this data is not enough, 
and finding the local structure of C12 around F12 would definitely require more work. 

In this paper we build on the work of Tadej and Zyczkowski [28j, by systematically 
studying the defect of complex Hadamard matrices. We will first recast the main results 
in [28] into our present algebraic geometric language. Then we will investigate the defect 
of the Fourier matrix Fq of an arbitrary finite abelian group G = Zjv^ x . . . x Z^v^- By 
using the same method as in [28], we will first obtain the following formula: 

\G\ 



d{FG) = Yl 



In order to compute this quantity, we use two observations, namely the fact that this 
quantity is multiplicative over the isotypic components of G, and the fact that for a 
product of p-groups, the number of elements of order < p^ is multiplicative as well. 
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These observations will lead to the following general formula, valid for any finite abelian 
group G, and which is our main result in this paper: 



diFc) = n 1 



Here the numbers oq 



(r-fc)afe_i + (ai + ...+afe_i)- 



k=l 

and ai < a2 < 



Gp = Zpai X ... X Zpor is the p-component of G, and [a]g 



depending on p, are such that 
= l + g + + + 

Finally, we will comment on the general question on whether the associated quantum 
algebraic objects, like the subfactor, planar algebra, or quantum permutation group, "see" 
or not the defect. Our remark here is that for a finite abelian group G we have: 

I 

k 



d{FG) = \im \im\y^ [ Xr 



Here the quantities Xr on the right are Diaconis-Shashahani type variables [H] asso- 
ciated to the embedding G C Sn, with = \G\. Now since for an arbitrary quantum 
permutation group G C these variables make sense, cf. [5], the above formula makes 
sense as well, and we will raise here the question regarding its possible vahdity. 

The paper is organized as follows: 1 is a preliminary section, in 2-3 we review the 
notion of defect, from the above-mentioned algebraic geometric point of view, in 4-5 we 
investigate the Fourier matrices, and 6 contains the probabilistic speculations. 

1. Complex Hadamard matrices 

We consider in this paper various square matrices over the complex numbers, H G 
M]\f{C). The indices of our matrices will usually range in the set {0, 1, . . . , — 1}. 

Definition 1.1. A complex Hadamard matrix is a matrix H G Mn{C) whose entries are 
on the unit circle, \Hij\ = 1, and whose rows are pairwise orthogonal. 

It follows from definitions that the columns of H are pairwise orthogonal as well. 
The main example is the Fourier matrix, based on the root of unity w = e^'^*/^: 





(I 


1 


1 


• 1 \ 




1 


w 






Fn = 
















^2(7V-1) 





There are many other complex Hadamard matrices, see [27] ■ One simple way of con- 
structing new examples is by taking tensor products, as for instance: 



F2®F2 



1 



1 



A 1 


1 1\ 


1 -1 


1 -1 


1 1 


-1 -1 


V -1 


-1 1/ 
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Observe that is nothing but the matrix of the discrete Fourier transform, over the 
cychc group Z^v. More generally, we have the following construction: 

Proposition 1.2. The Fourier matrix of a finite abelian group G — Zjvj x . . . x Z^y^ is 
the complex Hadamard matrix Fq = -F/Vi (8) . . . ® F^^ . 

Proof. For a product of groups G = G' x G" we have Fq = Fc ® Fq", and together with 
the above observation regarding F^, this gives the equality in the statement. □ 

As an example, for the Klein group Z2 x Z2 we obtain the above matrix F2® F2. 

Now back to the general case, observe that a matrix remains Hadamard when permuting 
rows and columns, or multiplying them by complex numbers of modulus 1. This is of 
course quite a trivial observation, and it is convenient to use the following notion: 

Definition 1.3. Two matrices are called '^equivalent" if one can pass from one to the other 
by permuting rows and columns, or multiplying them by complex numbers of modulus 1. 

Observe that each Hadamard matrix can be supposed, up to equivalence, to be in 
"dephased" form, in the sense that the first row and column consist of 1 entries only. 
Note that the Fourier matrices Fq = F^^ . . . (g) F^^ are by definition dephased. 

The Fourier matrices are of course quite well understood. The following general con- 
struction, inspired from [13], allows one, starting for instance from two Fourier matrices, 
to construct some quite non-trivial examples of complex Hadamard matrices: 

Definition 1.4. The deformation of a tensor product H ® K E Mnm{C), with matrix of 
parameters L G Mmxn{^), is H ®lK = {HijLajKab)ia,jb- 

Observe that for the "flat" matrix of parameters, Laj = 1, we obtain the usual tensor 
product H ® K. Observe also that in the above definition we can always assume L to be 
"dephased" , in the sense that its first row and column consist of 1 entries only. 

As a first example, by deforming the tensor product F2 ® F2, with the matrix of pa- 
rameters L = (j^ ^), we obtain the following complex Hadamard matrix: 



This matrix is in fact the only one at = 4, so let us take a closer look at it. As a 
first observation, at g G —1, —i} we obtain tensor products of Fourier matrices: 

Proposition 1.5. The matrices F22 o,re as follows: 

(1) Atq = lwe have Fl^ = -F2 F2. 

(2) Atq = -l we have Fl^ ^ F2 ® F2. 




/I 



1 



1 



1 \ 



(3) At q = ±i we have F22 
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Proof. The first assertion is clear, and the second one follows from it, by permuting the 
third and the fourth columns: 

/l 1 1 1 \ /l 1 1 1 \ 



-'^2,2 



-1 



1 



-1 



1 



V 



-1 



1 



-1 1 
-1 -1 

VI -1 1 -V 

As for the third assertion, this follows from the following computation 

/l 1 1 1 \ /l 1 1 1 \ 



2,2 



1 
1 



-1 ±i 
1 -1 
-1 



-1 



1 

1 



2,2 



Here we have interchanged the second column with the third one in the case g = i, and 
we have used a cyclic permutation of the last 3 columns in the case q = —i. □ 

As a second example, at = 6 we have two possible deformations of the Fourier matrix 
Fq = F2 ® = ^ F2. In dephased form, these matrices are given by: 



2,3 



(1 
1 

u 



1\ 

r 



'3 , 



3,2 



In terms of the parameter q 



-^2,3 



/I 

1 

1 
1 
1 

Vi 



3 3 

f 3 

1 1 

3 f 

f 3 



1 

r 
s 

-1 
— r 
—s 



1 

-1 



■ 2 
3 r 

js 

-1 

■ 2 
-3 r 



f9 

-^3,2 



b. j 




we have 






/I 


1 


1 


1 


1 


1 \ 


1 


-1 


r 


— r 


s 


—s 


1 


1 


J 


3 


3' 


3' 


1 


-1 


jr 


-jr 


fs 




1 


1 


3' 


3' 


3 


j 


VI 


-1 


■ 2 


■1 
-3 r 


js 


-js J 



There are many other examples of Hadamard matrices at A^ = 6, see [27]. Of particular 
interest are the following two matrices, due to Haagerup and Tao 



-"6 



/I 
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1 


1 


1 


1\ 


1 


-1 


i 


i 


—i 


—i 


1 


i 


-1 


—i 


q 


-q 


1 


i 


—i 


-1 


-q 


q 


1 


—i 


q 


-q 


i 


-1 


Vi 


—i 


-q 


q 


-1 


^ 1 



/I 

1 

1 
1 
1 

Vi 



1 

3 
1 



1 



333 
3^ 3^ 3 



33^33 
j^ j^ j 1 3 
j^ j 3^ 3 1/ 



The point with these matrices is that they are "regular", in the sense that any scalar 
product between distinct rows decomposes as a sum of quantities of type X]fc=i '^'^^1 with 
A G T and w = e^'^*/". Observe that the deformed Fourier matrices are regular too. 
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The complex Hadamard matrices of size < 5 were classified by Haagerup in [16] and 
the regular matrices at = 6 were classified in pj, and we have: 

Theorem 1.6. The complex Hadamard matrices of small order are as follows: 

(1) F2, -F3, F22, Fb '^'^e the only examples at N = 2,3,4,5. 

(2) F^2j Hg,TQ are the only regular examples at N = 6. 

Proof. The idea is that, the rows of our matrix being pairwise orthogonal, we must first 
understand the solutions of Oi + . . . + otv = 0, with |ai| = . . . = |aAr| = 1. 

At = 2,3,4 this is very simple: our equation must be, up to a permutation of the 
terms, a "trivial" equation of type a — a = 0, a + ja+ j'^a = Oora — a + 6 — 6 = 
respectively, with \a\ = \b\ = 1 and j = e^'^^l'^ . Now by using this observation, and trying 
to build a complex Hadamard matrix of size A^, this leads to F2, F3, only. 

At A^ = 5, 6 understanding the solutions of the above equation is a particularly difficult 
task. However, Haagerup was able to obtain the above result (1), see pJij. As for (2), this 
result is from j3jj, with the proof this time being also long, but purely combinatorial. □ 

2. Enveloping tangent spaces 

We study now the real algebraic manifold formed by all the N x N complex Hadamard 
matrices. Let us begin with some definitions: 

Definition 2.1. Let M^(C) C Mjv(C) be the set of matrices having 1 on the first row 
and column, and make Sn-i x •S'at-i act on it by permuting rows and columns. We set: 

(1) C]\i = M]\i{T) n \/NU]\f: the manifold of N x N complex Hadamard matrices. 

(2) = Cat n M^(C).- the submanifold consisting of dephased matrices. 

(3) Em = Dn/{Sn-i x Sn~i): the equivalence classes of matrices in Cn- 

Here, and in what follows, by "manifold" we will always mean "real algebraic mani- 
fold". Observe that Cn,Dn are indeed real algebraic manifolds, as being by definition 
intersections of such manifolds. As for Ej^, as defined above, this is just a set. 

Observe that we have surjective maps, as follows: 

Cn — > Dn En 

Here the first map is by definition obtained by "dephasing" the matrix, in the obvious 
way, and the second map is the canonical quotient map. Note that the dephasing map 
Cn ^ Dn is continuous, and is a covering of real algebraic manifolds. 

Let us first record a result coming from the classification results above: 

Proposition 2.2. The sets En with N small are as follows: 

(1) E2 = {F2}, Es = {F,}, E, = {F5}. 

(2) E^ = {Fl2\qeT}. 

(3) E, = {Fyq e T2} U {F^q e T'} U {Hl\q G T} U {T,} U X. 
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Proof. This is just a reformulation of Theorem 1.6 above, with X being by definition the 
set of equivalence classes of non-regular 6x6 complex Hadamard matrices. □ 

With this result in hand, we can now formulate a few basic observations: 

Proposition 2.3. The manifolds Cm,Dn in general not smooth, and not connected. 
Nor they are in general complex algebraic manifolds. 

Proof. All the assertions basically follow from Proposition 2.2, by using the canonical 
surjective maps ^ Djsi ^ Ep^ in order to get back to D]\[ and C^. 

Indeed, from \E2\ = 1 and from i?4 ~ T we obtain that C2, -D4 have real dimensions 3 
and 1, so they cannot have a complex algebraic manifold structure. 

Regarding now the smoothness claim, let us first look at D4. This manifold consists of 
|5'3 X S'al = 36 copies of G T} ~ T, that we will denote as follows: 



T(-''^) = {(7r,a)F|,|gG 



T} 



Here the element {n , a) E S3 x S3 acts as usual on M'^{C), with n permuting the rows 
labeled 1,2,3, and a permuting the columns labeled 1,2,3 (recall that, according to our 
usual conventions, the first row and column of our matrices are labeled 0). 

The point now is that some of these 36 copies of T intersect at g = 1, which prevents 
D4 from being smooth. For instance the copy of T associated to ((12), (12)) is: 



-^((12), (12)) 



/I 

1 



\1 -q 



1 \ 

-1 

-Q 
1J 



We can see from this formula that T*^*^^^^'*^^^^) is distinct from T*^*'^'*'^) = {-FI2I5' ^ "^l; 
but intersects it at g = 1, so is indeed not smooth. Similarly, C4 is not smooth either, 
because it consists of 36 copies of T x ~ T®, which intersect non-trivially at g = 1. 

For the connectedness question now, this follows from the fact that Tg and its conjugates 
are isolated in Dq, cf. [23], [28]. By lifting, Cq cannot be connected either. □ 

Regarding now the set En, it is not clear how to give it a structure of real algebraic 
manifold, or even a nice topological space structure. Observe that we have an embedding 
En C Mn(T) coming from the lexicographic order on Mn(T) induced by the [0, 2n) order 
on T. It is not clear whether this embedding produces a nice space or not. 

In order to get now some insight into the structure of Cn,Dn, we will compute the 
corresponding enveloping tangent spaces, coming from Definition 2.1 above: 

Definition 2.4. The "enveloping tangent spaces" ofCN,DN are defined by: 

(1) For H eCn we let ThCn = ^M^(T) n ThVNUn- 

(2) For H eDn we let ThDn = ThCn n ThM'^{C). 
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Here, and in what follows, we denote by ThM the tangent space to a smooth real 
manifold M, at a given point H G M. Observe that, Cn,Dn being not smooth manifolds 
in general, TnCjq^THDjq are of course not their "tangent spaces" in the usual sense. 

Tadej and Zyczkowski computed in [28] a certain vector space which coincides with the 
above enveloping tangent space. Their result, formulated in our terms, is as follows: 

Theorem 2.5. We have a canonical identification 

ThCn = |a G M^(M)| ^ H^kHjkiAk - Ajk) = 0, V«, j| 

and ThDn consists of the matrices A G T^Cn having on the first row and column. 

Proof. We use the notation Hij = Xij + iYy. We know that M]\f{T) is defined by the 
algebraic relations l-ffjjP = 1, with i,j G {0, 1, . . . , — 1}, and we have: 

Consider now an arbitrary vector v G ThM^IC), written as: 

'" = Y1 ^^^^^^ + f^^^"^^^ 

By taking the scalar product with the above quantity, we get: 

< V, d\Hij\^ >= 2{aijXij + /3ijY,j) 

Now since these scalar products vanish when we have aij = AijYij and Pij = —AijXij, 
for certain numbers Aij G M, we obtain the following formula: 



ThM^(T) = { J2Aij(Y,jX,,-XijY,, 



Aij G 



Let us compute now the subspace in the statement. We know that v NUn is defined by 
the algebraic relations < Hi, Hj >= N6ij, where Hq, . . . , H^^i are the rows of H. Since 
the relations < Hi, Hi >= N are automatic for matrices H G Mn{T), we can remove 
them from the computation. So, if for i 7^ j we let Lij = d < Hi, Hj >, then we have: 

fnCr, = {ve ThM^(T)| < V, L,, >= 0, ^ j} 
So, let us compute the scalar product appearing above. First, we have: 

Lij = d < Hi, Hj >= d i ^ HikHjk \ = ^ HikHjk + HjkHik 

\ k J k 

= + iyik){Xjk — iYjk) + i^jk — iyjk){^ik + iyik) 
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Consider now an arbitrary vector v G ThMn{T). According to the above formula of 
ThMi>4{T), we can write this vector in terms of a real matrix (A^fc), as follows: 

Ik 

Thus the above scalar products defining the space ThCn are given by: 

< V, Lij > = < AikiYikXik — XikYik) + AjkiYjkXjtt — Xj^Yj^), > 

k 

= ^ikiXik + iXik){Xjk + i^jk) + Ajk{Yjk — iXjjt){Xik — iYik) 

k 

= i ^ ^ AikHij^Hjk — AjkHjkHik 

k 

= i ^^(^jfc — Ajk)HikHjk 

k 

Thus we have reached to the description of Cat in the statement. The description of 
D^r is the statement follows as well, by intersecting with ThM'^{G). □ 

Let us look now at the dimensions of the spaces in Theorem 2.5: 

Definition 2.6. Associated to a complex Hadamard matrix H G Mjv(C) are: 

(1) The dephased defect d'{H) = dim{THDN). 

(2) The undephased defect d{H) = dim(TH-CAr). 

Observe that the dephased defect d'{H) is nothing but the quantity def{H) from [28] . 
simply called "defect" there. In what follows we will rather use the quantity d{H), which 
behaves better with respect to the various operations on complex Hadamard matrices. 

The dephased and undephased defect are related by the following formula: 

Proposition 2.7. d'{H) = d{H) -2N + 1. 

Proof. Consider the vector space M^(R) C Mn(M.) formed by the matrices of type Aij = 
tti + bj, with tti, bj G M. We claim that we have an isomorphism as follows: 

Indeed, the first remark is that for a matrix of type A^j = Oj + bj we have: 
HikHjk{Aik — Ajk) = ^ HikHjk{ai + bk — aj — bk) 

k k 

= {ai — aj) ^ HikHjk 

k 

= {ai - aj)6ij = 
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Thus we have M^(M) C ThCn- Also, we have M^(M) nTj^Z^jy = {0}, because tti + bj = 
if z = or j = imphes ai + bj = for any Finally, let A G ThCn, and decompose it 
as A = B + C, with Bij = Aio + Aoj-Aoo and Qj = Aij-Bij. Then we have B e M^(M), 

so C G ThCn, and since Cio = Coj = for any we conclude that C G ThDn- 

Summing up, we have proved the above vector space claim. Now since we have 
dim(M^(R)) = 2N — 1, we obtain the formula in the statement. □ 

The study of the defect is motivated by the following simple fact: 

Proposition 2.8. If d'{H) = then H is isolated inside D^. 

Proof. Indeed, if the matrix H E Dn was not isolated, then the tangent vector to any 
one-parameter family H'^ would belong to the space ThD^ = 0, contradiction. □ 

In general, it is not very clear what exact local information about Cat is encoded by 
the defect. For a detailed discussion here, we refer to the recent article [8]. 



We discuss in this section the various particular cases of Theorem 2.5. Let us first 
examine the product operations. We have here the following statement: 

Proposition 3.1. The enveloping tangent spaces are as follows: 

(1) For equivalent matrices H H' we have T c:^ T' . 

(2) For tensor products H = H' ® H" we have T' ® T" C T . 

Proof. These results, from |26], [28], follow as well from Theorem 2.5: 

(1) This is clear from definitions, because the permutations and the multiplication by 
scalars act as well, and in a compatible way, on the enveloping tangent spaces. 

(2) Indeed, for a matrix A = A' ® A" with A' G f' and A" G T", we have: 



3. Examples and comments 




fee 



fee 



II 



ac 



k c 



II 

be 



k c 



kc 




/ J ^ia^kc^ jh,kc-^jh,kc 
kc 



Thus we have indeed A G T, and we are done. 



□ 
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In terms of the defect, we obtain: 

Proposition 3.2. The undephased defect satisfies: 

(1) For equivalent matrices H H' we have d{H) = d{H'). 

(2) For tensor products H = H' H" we have d{H) > d{H')d{H"). 

Proof. These resuhs, once again from [2B], follow as well from Proposition 3.1. □ 

Observe that we don't have equality in the tensor product estimate, even in very simple 
cases. For instance if we consider two Fourier matrices F2, we obtain: 



Here the number 10 comes from the general formula of (i(-Fb) explained in section 5 
below, in the particular case G = Z2 x Z2, or from Proposition 3.5 below, at g = 1. 

The problem of finding upper bounds for the defect of a tensor product was investigated 
in detail by Tadej in [26]. We would like to raise here the following related question: 

Problem 3.3. How does the defect of a deformed tensor product, d{H ®l K), vary with 
the matrix of parameters L G MmxnO^) ? 

This problem looks quite complicated, even at small values of M, A^. As a first obser- 
vation, the equations of the enveloping tangent space are: 



There is no obvious trick that can be applied here. Note that there is no reason for a 
"transport formula" of type d{H®LK) = d{H®K) to hold, even in simple cases. Indeed, 
the L-deformation procedure "deforms well" the matrix H ^ K, but maybe not the other 
NM X NM complex Hadamard matrices, which might happen to be around. 

However, there are several special cases where the problem can be solved, for instance 
by combining the formulae for Fourier matrices with the following observation: 

Proposition 3.4. We have F^m — ®l Fm, where Laj = w"-^ , with w = e^'^*/^^'^. 

Proof. Indeed, by using w^'^ = 1, we obtain the following formula: 



d{F2 ® F2) = 10 > 9 = d{F2Y 




kc 



{Fn ®L FM)ia,jb — {FN)ijLaj{FM)ab 



^Mij+aj+Nab 
XNi+a){j+Nb) 



Here F; 



— iF^;^f)ia,jb 

is a certain matrix which is equivalent to F^m, and we are done. 



□ 



Let us discuss now the simplest case of the problem, N = M = 2. We will work out 
everything in detail, as an illustration for how the equations in Theorem 2.5 work. 



12 TEODOR BANICA 

Proposition 3.5. We have d{F2 2) = 10 at g = ±1, and d{F2 2) = 8 at q ^ ±1. 
Proof. Our starting point are the equations in Theorem 2.5, namely: 

HikHjk{Aik — Ajk) — 

h 

Since the i > j equations are equivalent to the i < j ones, and the i = j equations are 
trivial, we just have to write down the equations corresponding to indices i < j. And, 
with ij — 01, 02, 03, 12, 13, 23, these equations are: 



{Aoo- 


Aio) - {Aoi - 


An) + q{Ao2 - 


A12) - q{Ao3 


-A13) 


= 


{Aoo 


- A2q) + {Aoi 


-A21)- {A02 - 


- A22) - {Aq3 


-A23) 


= 


{Aoo- 


A30) - {Aoi - 


A31) - q{Ao2 - 


A32) + q{Ao3 


-^33) 


= 


{Aio- 


A20) - {All - 


A21) - q{Ai2 - 


A22) + q{Ai3 


-^23) 


= 


{Aio 


- Aso) + {An 


-Asi)- {A12 - 


- A^2) - {Ai3 


-^33) 


= 


{A20- 


A30) - {A21 - 


A31) + q{A22 - 


A32) - q{A23 


-^33) 


= 



Assume first q ±1. Then q is not real, and appears in 4 of the above equations. But 
these 4 equations can be written in the following way: 



(^0 


-Aoi) 


-{Aio 


- All) + q{{Ao2 


-A03) 


-{A12 


-Mz)) = 





(^00 


-Aoi) 




- Asi) - q{{Ao2 


-A03) 


-{A32 


-A33)) = 





(A-io 


-All) 


- (^20 


- A21) - q{{Ai2 


-A13) 


- {A22 


-A23)) = 





{A20 


-A21) 


- (^30 


-A3i)+q{{A22 


-A23) 


- {A32 


-A33)) = 






Now since the unknowns are real, and q is not, we conclude that the terms between 
braces in the left part must be all equal, and that the same must happen at right: 

^00 ~ ^01 — Aio — All = A20 — A21 = A30 — A31 
A02 — Ao3 = A12 — Ai3 = A22 — A23 = A32 — A33 
Thus, the equations involving q tell us that A must be of the following form: 

/ a a + X e + y e\ 

b b + x f + y f 

c c+x g+y g 

yd d + X h + y hj 

Let us plug now these values in the remaining 2 equations. We obtain: 

a — c + a + X — c — X — e — y + g + y — e + g = 
h-d + h + x-d-x- f -y + h + y- f + h = 

Thus we must have a + g = c + e and b + h = d + f , which are independent conditions. 
We conclude that the dimension of the space of solutions is 10 — 2 = 8, as claimed. 
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Assume now q = ±1. For simplicity we set g = 1, by using Proposition 3.1, and we use 
as well Proposition 2.7, which tells us that it is enough to compute the dephased defect. 
The dephased equations, obtained by setting Aio = Aqj = in our system, are: 

^11-^12 + ^13 = 

-A21 + A22 + A23 = 

^31+^32-^33 = 

-An + A21- A12 + A22 + Au- A23 = 
An - Asi - Au + A32 - An + A33 = 
-A21 + A31 + A22 - A32 - A23 + A33 = 

The first three equations tell us that our matrix must be of the following form: 

(a a + b b \ 
c+d c d \ 
e f e + fj 

Now by plugging these values in the last three equations, these become: 

—a + c + d — a — b + c + b — d = 
a — e — a — b + f — b + e + f = 
-c-d+e+c-f~d+e+f = 

Thus we must have a = c, b = f,d = e, and since these conditions are independent, 
the dephased defect is 3, and so the undephased defect is 3 + 7 = 10, as claimed. □ 

Let us discuss now the reformulation of Theorem 2.5, for certain special classes of 
complex Hadamard matrices. The simplest situation is that when we have a usual 
Hadamard matrix, H G M7v(±l). The combinatorics of H is encoded into the array 
e G M7VxArx7v(±l) given by eijk = HikHj^, that we will call "design" of H. See pQ, [T5] . 

Proposition 3.6. In the real case H G M7v(±l) the system of equations for T is 

^ ] ^ijkj^ik — ^jk) = 
k 

where e G MArxArxAf(±l); given by Sijk = HikHjk, is the "design array" of H. 

Proof. This is just an observation, clear from definitions. □ 

It is not clear what the spectral properties of e are, and how the above formulation 
can be improved. Let us also mention that of particular interest would be an extension 
of the above observation to the root of unity case [T3]. However, as shown by Lam and 
Leung in [22], the "design" of a Butson matrix can be something of a very complicated 
nature, unless the "regularity conjecture" in [3] holds indeed. One way of overcoming 
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these problems would be by looking directly at the matrices which are regular in the 
sense of pj. But the corresponding "design matrices" are not axiomatized yet. 

Let us discuss as well the circulant case. Here our matrices must be "Fourier-diagonal" , 
of type H = \/NFQF*, with F = F]\[/\/N and with Q diagonal over T. The problem of 
finding the diagonal vectors Q G having the property that H = \fNFQF* is indeed 
Hadamard is a subtle Fourier analysis problem, cf . [7] , [TT] , [12] , [IZ] . 

In Fourier notation, our system of equations is as follows: 

Proposition 3.7. In the circulant case the equations for T are 

klr 

where Q E is the eigenvalue vector of U = H/ \fN . 
Proof. First, with if^ = Cj_i, the system of equations for T becomes: 

''^^^Ck-iCk-jiAik — Ajfc) = 

k 

Now since we have C = FQ, where Q G is the eigenvalue vector, we get: 

Mr 

By rearranging the terms, this gives the formula in the statement. □ 

Once again, there is no obvious trick that can be applied to the equations. The circulant 
problem seems to require a delicate case-by-case analysis, and we have no results. 

We would like to end this section with a few theoretical observations, regarding the 
general case. First, we have the following linear algebra interpretation of the defect: 

Proposition 3.8. The undephased defect d{H) is the corank of the matrix 

{Re{HibHjb) i < J 
Im{HihHjb) i > j 
* i = j 

where * can be any quantity (its coefficient being anyway). 

Proof. The matrix of the system defining the enveloping tangent space is: 

Xij,ab = i^ia — SjajHibHjb 

However, since we are only looking for real solutions A G Mjv(M), we have to take into 
account the real and imaginary parts. But this is not a problem, because the (ij) equation 
coincides with the (ji) one, that we can cut. More precisely, if we set Y as above, then we 
obtain precisely the original system. Thus the defect of H is indeed the corank of Y. □ 
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As an illustration, for the Fourier matrix Fn we have the following formula, where 
e{i,j) e {—1, 0, 1} is negative ii i < j, null for i — j, and positive for i > j: 

Observe in particular that for the Fourier matrix F2 we have: 

/ \ 

-111-1 
\ / 

Here the corank is 3, but, unfortunately, this cannot be seen on the characteristic 
polynomial, which is i-*(A) = A^. The problem is that our matrix, and more precisely its 
middle 2x2 block, is not diagonahzablc. This phenomenon seems to hold in general. 

Finally, we have as well the following key question: 

Question 3.9. Does the profile matrix of H, namely 

k 

determine the enveloping tangent space ThCn? 

The point here is that the profile matrix M e Mn2 (C) is the one which produces the 
subfactor, planar algebra, or quantum permutation group associated to H E Mjv(T). 
Thus, the above question is of great importance for understanding the relation between 
the "geometric" and "quantum" invariants of the complex Hadamard matrices. 

An idea for parametrizing the system in terms of the profile matrix M might come 
from the solution for the Fourier matrix, explained in section 4 below. Indeed, if we write 
P — AH/N, so that A — PH*, then our equations become: 



HksHikHjk ^(Pjs — Pjs) — 

k s 



Observe the similarity between the matrix on the top left and M. However, this does 
not solve our problem. We will be back to our problem in section 6 below. 

4. Abstract Fourier matrices 

We recall that the tensor products of Fourier matrices are exactly the Fourier matrices 
of the finite abelian groups. More precisely, with G = Zjvi x . . . x Zjv^, we have: 



Fa — -Fj 
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The defect of -Fjv was computed by Tadej and Zyczkowski in [2H] . The first part of their 
proof works in fact for abstract Fourier matrices as well, and we have: 

Proposition 4.1. For a Fourier matrix F = Fq, the matrices A G TfCn, with N = \G\, 
are those of the form A = PF* , with P G Mn{C) satisfying 

where the indices i,j are by definition taken in the group G. 
Proof. We use the system of equations in Theorem 2.5, namely: 

FikFjk{Aik — Ajk) = 

k 

Now assume F = F/Vi ® ■ ■ ■ ® F^^^ , so that with Wk = e^'^'^l^ we have: 

i^n...W....> = (^^J^^^^ • • 
With = A^i . . . and w = e^'^^l^ , we obtain: 

Thus the matrix of our system is given by: 

- f ('l-jl)fcl I (ir-jr)kr \ »r 

F- ■ 1 1 F- ■ u 1 = vA ^1 J 
Now by plugging in a multi-indexed matrix A, our system becomes: 

Now observe that in the above formula we have in fact two matrix multipliations, so 
our system can be simply written as: 

(^-P^)n...ir,n-ji---ir-ir ~ (^-^)ii---jr,n-ji-.-v-jr ~ 

Now recall that our indices always have a "cyclic" meaning, so they belong in fact to 
the group G. So, with P = AF, and by using multi-indices, our system is simply: 

p. . = P 

With i = I + J, j = I we obtain the condition Pi+j,j = Pu in the statement. 
In addition, A = PF* must be a real matrix. But, if we set Pij = Pi-j-, we have: 

iP'-P*)il...ir,jl...jr ~ ^ ] Pii...ir,ki---kr-Fjl—jr,kl---kr 

= ^ ] Pil...ir,-k\...-kr{,F )-kl...-kr,jl...jr 
kl...kr 

— {PF )ii...ir,jl...jr 
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Thus we have PF* — PF*, so the fact that the matrix PF* is real, which means by 
definition that we have PF* — PF*, can be reformulated as PF* — PF*, and hence as 
P — P. So, we obtain the conditions Pij — Pi-j in the statement. □ 

Theorem 4.2. The undephased defect of an abstract Fourier matrix Fq is given by: 

Proof. According to Proposition 4.1 above, the undephased defect d{Fc) is the dimension 
of the real vector space formed by the matrices P e M]sf{C) satisfying: 

Here, and in what follows, the various indices . . . will be taken in G. Now the point 
is that, in terms of the columns of our matrix P, the above conditions are: 

(1) The entries of the j-th column of P, say C, must satisfy Cj = Cj+j. 

(2) The (— j)-th column of P must be conjugate to the j-th column of P. 

Thus, in order to count the above matrices P, we can basically fill the columns one by 
one, by taking into account the above conditions. In order to do so, consider the subgroup 
G2 — {3 e G|2j = 0}, and then write G as a disjoint union, as follows: 

G = GaUXU (-X) 

With this notation, the algorithm is as follows. First, for any j e G2 wc must fill the 
j-th column of P with real numbers, according to the periodicity rule Q = Q+j. Then, 
for any j e X we must fill the j-th column of P with complex numbers, according to the 
same periodicity rule Cj = Cj+j. And finally, once this is done, for any j G X we just 
have to set the (— j)-th column of P to be the conjugate of the j-th column. 

So, let us compute the number of choices for filling these columns. Our claim is that, 
when uniformly distributing the choices for the j-th and (— j)-th columns, for j ^ G2, 
there are exactly [G :< j >] choices for the j-th column, for any j. Indeed: 

(1) For the j-th column with j G G2 we must simply pick real numbers subject to 
the condition Q = Cj+j for any i, so we have indeed [G :< j >] such choices. 

(2) For filling the j-th and (— j)-th column, with j ^ G2, we must pick N complex 
numbers subject to the condition = Cj+j for any i. Now since there are [G :< j >] 
choices for these numbers, so a total of 2[G' :< j >] choices for their real and imaginary 
parts, on average over j, — j we have [G :< j >] choices, and wc arc done again. 

Summarizing, the dimension of the vector space formed by the matrices P, which is 
equal to the number of choices for the real and imaginary parts of the entries of P, is: 

diFc) = J2[G :< J >] 
jeG 

But this number is exactly the one in the statement, and this finishes the proof. □ 
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The above results suggest a number of concrete and abstract computations, to be done 
in the next sections. For the moment, let us just record the following conjecture: 

Conjecture 4.3. For a regular Hadamard matrix H G MAr(C) we have 

ThCn = C ■ [ThCn]q 

where [ThCn]q consists of the matrices A G ThCn having rational entries. 

More precisely, this result seems to hold for the Fourier matrices, and we conjecture 
that this result holds in fact in general, under the "regularity" assumption in [3]. 

For usual Fourier matrices F/v, the result definitely holds aX N = p prime, because the 
minimal polynomial of w over Q is simply P = 1 + w + . . . + w^~^. The case N = has 
also a simple solution, coming from the fact that all the px p blocks of our matrix A can 
be shown to coincide. In general, this method should probably work for N = p^, or even 
for any G N. However, since we don't have a complete proof here, the Fourier matrix 
statement should be regarded as being part of the above conjecture. 

5. Finite group calculations 

In this section we complete the computation of the defect of Fq, by using the formula 
in Theorem 4.2. It is convenient to consider the quantity 6{Fg) = \G\^^d{FG), which 
behaves better, and to try to compute it in the abstract group framework: 

Definition 5.1. Associated to a finite group G is the following quantity: 

As a first example, consider a cyclic group G = Zj^i, with N = p"" power of a prime. 
The count here is very simple, over sets of elements having a given order: 

5(Zp.) = 1 + (p - l)p'^ + {p" - p)p-^ + . . . + (p'^ - p''-^)p'^ = 1 + a - - 

p 

In order to extend this kind of count to the general abelian case, we use two ingredients. 
First is the following lemma, which splits the computation over isotypic components: 

Lemma 5.2. For any finite groups G, H we have: 

6{G xH)> 6{G)6{H) 

In addition, if {\G\,\H\) = 1, we have equality. 

Proof. Indeed, we have the following estimate: 



ord{g,h) ^ [ord{g),ord{h)] 



^-^ ord[g) ■ ord[n) 



gh 
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Now in the case (|G|,|i?|) = 1, the least common multiple appearing on the right 
becomes a product, [ord{g) , ord{h)] — ord{g) ■ ord{h), so we have equality. □ 

Proposition 5.3. For a finite abelian group G we have 

5{G) = Y[5{G,) 

V 

where Gp with G — XpGp are the isotypic components of G. 

Proof. This is clear from Lemma 5.2, because the order of Gp is a power of p. □ 

The second ingredient concerns the p-groups, and is as follows: 
Lemma 5.4. For the p-groups, the quantities 

Ck = G G\ord{g) < 
are multiplicative, in the sense that Ck{G x H) = Ck{G)ck{H). 
Proof. Indeed, for a product of p-groups we have: 

Ck{GxH) = ^{{g,h)\ord{g,h)<p''} 

= #{(g,h)\ord(g)<p',ord(h)<p'} 
= M9\ord{g) < p'U{h\ord{h) < p'} 
We recognize at right Ck{G)ck{H), and we are done. □ 
Let us compute now 5 in the general isotypic case: 
Proposition 5.5. For G = Zp^i x . . . x Z^a^ with ai < a2 < ■ ■ ■ < ar we have 

r 

5{G) = 1 + Y^pir-k)a,-^+{a.+-+a,-^)-l^pr-k+l _ ^^^^^ _ 
k=l 

with the convention ao = 0, and by using the notation [a]q = l + q + q'^ + ... + 
Proof. First, in terms of the numbers c^, we have: 



S{G)^1 + J2 



Ck — Ck-l 



k>l ^ 

In the case of a cyclic group G = Zpa we have c^ — p™™^^'"). Thus, in the general 
isotypic case G — l^pix x . . . x Zpar we have: 

_ ^min(A;,ai) miii(fc,ar) _ ^min(fc,ai)+...+min(fc,ar) 

L^k — P . . . p — p 

Now observe that the exponent on the right is a piecewise linear function of k. More 
precisely, by assuming ai < a2 < . . . < Or as in the statement, the exponent is linear on 
each of the intervals [0, ai], [ai, 02], . . . , [otr-i, o,r]- So, the quantity S{G) to be computed 
will be 1 plus the sum of 2r geometric progressions, 2 for each interval. 
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In practice now, the numbers Ck are as follows: 

Co = 1, Ci = p'', C2 = /^ . . . , Cai = P™' , 

^ _ „ai+(r-l)(ai+l) _ ai + (r-l)(ai+2) _ ai + (r-l)a2 

•--oi + l F 1 F 1 • • • 1 '-a2 F i 

„ — „ai+02 + (r--2)(02+l) „ _ „ai+02+(r--2)(o2+2) _ ^oi+a2+(r-2)a3 

(-02+1 — ) 1-02+2 — F 1 ■ ■ ■ 1 i-aa — F i 



„ _ „ai+...+ar-i+(ar-i + l) „ _ „ai+...+a,._i + (ar_i+2) _ ai + ...+ar 

^ar-i+1 — ! — F 1 ■ ■ ■ 1 — F 

Now by separating the positive and negative terms in the above formula of 5{G), we 
have indeed 2r geometric progressions to be summed, as follows: 

-{p-^ + + p^""^ + . . . + 

+ _^p(''-l)(»i+2)-2 _^ _ _ _^^ai+(r-2)a2^ 

_|_p(»'-l)(ai+l)-2 _|_ _ _ _ _|_ pai+(r-l)(a2-l)-a2'j 

+ . . . 

_j_^pai+---+ar-i _|_ pO-i + .-.+ar-i _|_ _j_ pai+---+ar-l ^ 
_^pai+---+ar-l-l _j_ pQl+.-.+ar-l-l _j_ _j_ pOi+.-.+a^-i-l j 

Now by performing all the sums, we obtain: 

(r-l)ai _ 1 

5{G) = l+p-\f-l) \,_,_^ 

„(r-2)(a2-ai) _ i 

, (r-2)ai+(ai-l)/ r-1 _ -^\P t 

^ V/^ > pr-2_l 

.^pC?— 3)a2+(ai+a2-l)^^r— 2 _ -j^-j 



p(r-3)(a3-a2) _ ]^ 



pr-3 _ 1 
+ ... 

+p"^+-+"'-^"'(p-l)(a.-a,_i) 

By looking now at the general term, we get the formula in the statement. □ 

Let us go back now to the formula in Theorem 4.2. By putting it together with the 
various results in this section, we obtain our main result in this paper, namely: 

Theorem 5.6. For a finite abelian group G, decomposed as G = XpGp, we have 
d{FG) = \G\ n ( 1 + - l)[afe - ak-i]pr-, 

p \ k=l 

where Qq — and Oi < 02 < . . . < Or T depending onp) are such that Gp — Zpai x . . . x Zpa,. . 
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Proof. Indeed, we know from Theorem 4.2 that we have d{FG) = \G\6{G), and the result 
follows from Proposition 5.3 and Proposition 5.5. □ 

As a first illustration, we can recover in this way the defect computation in [28] : 

Corollary 5.7. The undephased defect of a usual Fourier matrix Fn is given by 

d{FN)=Nf[(^l + ai-^^ 
where N = . . . p'^" is the decomposition of N into prime factors. 

Proof. The underlying group here is the cyclic group G = Z^, whose isotypic components 
are the cyclic groups Gp^ = Zp^, . By applying now Theorem 5.6, we get: 

s 

d{Fr,) = Nl[{l+pr\p,-l)a,) 

i=l 

But this is exactly the formula in the statement. □ 
As a second illustration, for the group G = Zp^i x Zpa2 with ai < a2 we obtain: 
diFc) = p''^+'='%l+p-\p^-l)[a,]p + p'^^-\p-l)ia,-a,)) 

p —\ 

In particular at j9 = 2 and ai = 02 = 1 we obtain that the defect of the Fourier matrix 
Fg = -F2 ® F2, already known from Proposition 3.5 to be 10, is indeed: 

diFc) = 2(2 + 3 + 0) = 10 

Finally, let us mention that for general non-abelian groups, there doesn't seem to be 
any reasonable algebraic formula for the quantity 6{G). As an example, consider the 
dihedral group D^, consisting of symmetries and A^ rotations. We have: 

6{Dm) = y + 6{Zn) 

Now by remembering the formula for Zn, namely 5{Zn) = J|(l +p~'^{pi — l)cii), it is 
quite clear that the A^/2 factor couldn't be incorporated in any nice way. 

6. Probabilistic speculations 

We have seen in the previous section that the defect of a Fourier matrix Fq can be 
computed by doing some explicit calculations in the associated group G. In this section 
we speculate on a possible extension of this method, by using quantum groups. 

The story here goes back to Popa's paper |2l], who discovered that the "orthogonal 
MAS A" condition A _L UAU* holds inside the N x N matrices precisely when H = 
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\fNU is Hadamard. Such orthogonal MASA are known to produce subfactors, and Jones 
subsequently came with a complete study of the problem, from the point of view of 
statistical mechanics [IH], general subfactor theory [2T], and planar algebras [20] . 

Thanks to some general Tannakian correspondences, it was realized in the late 90's that 
these subfactors come in fact from certain quantum groups. More precisely, each complex 
Hadamard matrix H e MAr(C) produces a certain quantum permutation group G C 5^, 
and the subfactor associated to H can be understood in terms of G. 

The construction H G can be summarized as follows: 

Definition 6.1. Let H G Mn{G) he a complex Hadamard matrix. 

(1) Set = Hi/ Hj, where ifi, . . . , Hj^ G are the rows of H . 

(2) Let Pij G MAr(C) be the orthogonal rank one projection on 

(3) Define a representation tt : C{S^) — )■ MAr(C) by T^iuij) = Pij. 

(4) Consider the factorizations of type n : C{S^) — t- C{G) — t- Mjv(C). 

(5) Let G C be the minimal object producing such a factorization. 

In this definition is the quantum permutation group of Wang [3lJ, known to be a 
compact quantum group in the sense of Woronowicz [32]. For technical details regarding 
this construction we refer to [5], [0], and for a recent survey on the subject, to [2]. 

The basic result here, which is for interest for us, is: 

Proposition 6.2. The quantum permutation group associated to the Fourier matrix Fa 
of a finite abelian group G is nothing but the group G itself, acting on itself. 

Proof. The idea here is that examining the definition of vr leads to an obvious factorization 
of type TT : G{S^) G{G) Mn{C), with N = \G\. Now since G is classical, the 
minimality property of this factorization is not hard to establish. See [3]. □ 

This result, when combined with those in the previous section, raises the idea of com- 
puting the defect of an arbitrary complex Hadamard matrix H G Mjv(C) by using the 
associated quantum permutation group G C S'^. We do not know if this is possible: 

Question 6.3. Does the quantum permutation group G C see the defect of the complex 
Hadamard matrix H G M]\f{C) ? And if so, how exactly? 

In what follows we will just speculate on the second question. We recall from section 4 
that in the case of Fourier matrices we have indeed a formula, namely: 

The problem is that the quantity on the right is not available for arbitrary quantum 
groups G C S"^, simply because these quantum groups are abstract objects, having no 
elements g & G. So, let us pause now from our defect investigation, and remember what 
exact numeric quantities are available, for an arbitrary quantum group G C Sj^. 
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The answer comes from Woronowicz's paper [32], who developed there an analogue of 
the Peter- Weyl theory for the compact quantum Lie groups. In particular, Woronowicz 
established an analogue of the following key representation theory formula: 



Now in the case of the usual symmetric group Sn C On, the character Tr : G — > M 
on the left is nothing but the number of fixed points. So, as a conclusion, we can say 
that "a quantum permutation group G C Sjj- doesn't exist as a concrete object, but the 
fixed point statistics Tr : G — )■ M does exist as a noncommutative random variable, and 
its moments can be computed by using representation theory methods" . 

As an example, a well-known result in classical probability tells us that in the limit 
iV — )■ oo, the character Tr : N follows the Poisson law. By using representation 

theory methods, one can prove that in the limit N ^ oo the character Tr : — ^ M 
follows a free Poisson law, in the sense of free probability theory [BUj. See pj. 

In view of this observation, the following problem appears: 

Problem 6.4. In the case of Fourier matrices H = Fq, can one recover the defect d{H) 
from the fixed point statistics Tr : G — )■ N o/ the underlying group ? 

The answer here is of course "yes" , due for instance to the following formula: 

ord{g) = min{r G N|Tr(5f'~) = A^} 

In order to deduce a formula for the defect, it is convenient to use the normalized trace 
tr = Tr/N. The first remark is that the "min" can be replaced by two parameters: 



Indeed, with k ^ oo only the maximal values of tr, namely tr = 1, contribute to the 
computation of the limit, and these appear with a frequency of l/ord{g). 
Now by summing over all the group elements, we obtain: 



The point now is that the variables Xrid) = tr{g^) are familiar objects, called Diaconis- 
Shashahani type variables for G. See [H]. Thus, we obtain the following result: 







r=l geG 
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Theorem 6.5. For G C Sn abelian acting on itself, with N = \G\, we have 

d(FG) = iV2lim limy V / Xr 

r=l ''^ 

where Xr{g) = tr{g^) are the Diaconis-Shashahani type variables for G. 

Proof. This is clear from the above considerations. □ 

The point now is that for an arbitrary quantum permutation group G C the vari- 
ables Xr are also available, cf. [5]. So, we have the following question: 

Question 6.6. Given a complex Hadamard matrix H G M7v(C), is there a formula for 
d{H) in terms of the variables Xr over the associated quantum group G d ? 

Finally, let us mention that the planar algebra associated to H is known to be given 
by Pk = Fix{u®^), so that the planar algebra dimensions are: 

dimPfe = Xi 
Jg 

As a conclusion, it might happen that the planar algebra dimensions don't see the 
defect, but the Diaconis-Shashahani type variables, taken together, do see it. 
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